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ABSTRACT 

We investigate a relativistic fluid jet driven by radiation from a shocked accre¬ 
tion disc around a non-rotating black hole approximated by Paczyhski-Wiita 
potential. The sub-Keplerian and Keplerian accretion rates control the shock 
location and therefore, the radiation field around the accretion disc. We com¬ 
pute the radiative moments with full special relativistic transformation. The 
effect of a fraction of radiation absorbed by the black hole has been approxi¬ 
mated, over and above the special relativistic transformations. We show that 
the radiative moments around a super massive black hole are different com¬ 
pared to that around a stellar mass black hole. We show that the terminal 
speed of jets increases with the mass accretion rates,synchrotron emission of 
the accretion disc and reduction of proton fraction of the flow composition. To 
obtain relativistic terminal velocities of jets, both thermal and radiative driv¬ 
ing are important. We show for very high accretion rates and pair dominated 
flow, jets around super massive black holes are truly ultra-relativistic, while 
for jets around stellar mass black holes, terminal Lorentz factor of about 10 
is achievable. 

Key words: Black Holes, Jets and outflows, Hydrodynamics, Radiation 
dynamics, Shock waves 


1 INTRODUCTION 

Astrophysical jets are ubiquitous, as they are associated with many classes of astrophysical 
objects such as active galactic nuclei (AGN e.g., M87), young stellar objects (YSO e.g., 
HH 30, HH 34), X-ray binaries {e.g., SS433, Cyg X-3, GRS 1915-1-105, GRO 1655-40) etc. 
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However, only jets around X -ray binaries 


and AGN like 3C273, 3C345 f Zensus et al. 


ike GRS1915+105 


Mirabel fc Rodrigueall994j) 


1995), M87 (Biretta 1993) etc are relativistic. In 


this paper we concentrate on relativistic jets. Since it is conjectured that stellar mass black 
hole reside at the heart of microquasars and those of the super-massive variety dictates the 
dynamics of the AGNs/quasars, therefore it implies that a jet has to originate from the 
accreting matter itself, since black holes (hereafter, BH) have neither hard surface nor any 
atmosphere. Interestingly, simultaneous radio and X-ray observations of microquasars show a 
very st rong correlation between the spectral s t ates of the accretion disc and the associated jet 


states (IGallo et. al. 


2003 


Fender et al. 


2010 : 


Rushton et al. 


2010), which reaffirms the fact 


that jets do originate from the accretion disc. In addition, recent observations have shown 
that jets originate from a region which is less than 100 Schwarzschil d radii (rg) around the 

2 OI 2 I) . which imply that the 


unresolved central object fjjunor et. al. 


1999 


Doeleman et. al. 


entire disc may not participate in production of jets, but only the central region of the disc 
is responsible. 

Since jets are supposed to originate very close to the central object, the plasma at the 
base should be hot and is expected to be fully ionized. This hot outflowing plasma will also be 
in the intense radiation field from the accret ion disc. A number of scientists have studied the 

G98nl l g _ g d partic 


interaction of radiation with jets. 


Icke 


and gas flow in the radia tion field of an underlying Keplerian disc fjShaknra fc Snnyaevlll973l ) 


Sikora fc WilsonI fll98ll) studied the interaction of the particle jets with t he radiation field. 


in the funnel like region of a thick accretion disc flPaczvhski fc Wiitalll980[l . Ignoring gravity 
and for normal electron-proton or e~ — p'^ plasma jets, the authors achieved a terminal speed 
of around vt ~ 0.4c, and obtained ter mina l Lorentz factor y'r ~ 3 for electron-positron or, 
e~ — e"*" plasma. In a seminal paper, lickel fjl989l ) showed that for particle jets above an 
infinite Keplerian disc, the radiation drag ensures an upper limit of terminal speed, which 
the author termed it as ‘magic speed’ and which turned out to be Umagic = 0.45c. The much 
‘vaunted’ magic speed is actually the so-called equilibrium speed of jet plowing through 
a radiation field near the Keplerian disc surface. It may be noted that, equilibrium speed 
(ueq) is the speed of the jet at which the radiation force becomes zero and for speeds above 
which, radiative deceleration sets in. Equilibrium speed arises due to the presence of the 
radiation drag. Radiation drag is significant for a radiation field due to an extended source. 
On the other hand, terminal speed {vt) is the speed of the jet at which the total force on 
the jet approaches zero. Therefore, the jet may achieve vt only at large distances away from 
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Radiatively and thermally driven jets 3 
the central object, wh ile it may reach Ueq at distances much closer to the central object. 


Melia fc Konig. 


I (l3) 


considered jets which are accelerated to ultra-relativistic speed way 


above the local Ueq. This causes the radiation drag t o become effective and decelerate the jet 

(119961 1 concluded that for the disc they chose, 


to terminal velocity vt ~ 0.995. 


Sikora et al. 


the maximum possible of Lorentz factor is ^ 4. The Japanese group led by Jun Fukue made 
very important contribution to the research of tr ansonic outflows, both in the relativistic, as 
well as, in non-relativistic domain. iFukuel (1199611 extended Icke’s work, and studied particle 
jets away from the axis, although like Icke, considered near disc approximation for the 
radiation field above it. The actual terminal jet velocity achieved was vt ~ 'Wmagio which is 


expected. However, the main prob 
jet and thereby spreading the jet. 


em is t 


Pukue 


re ten dency of the radiation field to spin up the 
19991) then s tudied jets con f ined 


in order to arrest the spreading and collimate the jet. 


ry disc corona 


Hirai fc Fukud (120011) on the other 


hand, computed the radiation field due to a Keplerian disc governed by Newtonian gravity, 
Schwarzschild gravity and Kerr gravity. The strength of the radiation held above the disc, 
described by Schwarzschild gravity is 50% lesser than that due to Newtonian gravity. But the 
strength of the radiation held above a disc governed by a Kerr type gravity is 10 times higher 
than that due to the Newtonian gravity, maki ng radiatively drive n jets easier to blow for a 


rotating black hole. In a very interesting paper. 


Fukue et al 


(I 2 OOII) considered a hybrid disc. 


consisting of outer K eplerian disc and inner advection dominated accretion how or ADAF 
(INarayan et a/.lll997l ) type accretion solutions. Since ADAF is dimmer, so the inner region 
from which the electron-positron jet is assumed to emerge in this paper, does not contribute 
in the radiation held. Such a scenario do produce jets with terminal Lorentz factor 'yx ~ 2, 
and the radiation held from the outer Keplerian disc also helps in collimation. 

Along with the various disc moc 
on ad vective discs was initiated by 


els like Keplerian disc, thick disc, A 

DAF, investigations 

Liang & Thompson 

(l9aQ 

); 

Fukue 

(19§7 

); 

Chakrabarti 


(119891 ). Such a disc can admit smooth solutions, and for the right choice of parameters it 
may harbour shock transition. This disc model was extended by considering injection of a 
mixture of matter with Keplerian angular momentum and sub-Keplerian angular momen¬ 
tum. The portion of disc which is termed as sub-Keplerian disc (SKD), sandwiches the 
Keplerian disc (KD) from the top and the bottom (see Fig. 1). SKD may undergo shock 
transition and at the shock, KD termin ates due to extra heating in the post-shock disc 
or PSD (jChakrabarti fc TitarchnkI Il995l ). Although it was initially proposed as an elab¬ 
orate and contrived model solution, but was recently conhrmed by numerical simulation 
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4 Vyas et al. 


studied 

ing jet ( 

by 

Ghattooadhvav & Ghakrabarti 

(2000a 

b, 

2002lj 

Ghattooadhvav et al. 

2004; 

Ghattooadhvav 

20051 


The investigation of outflow- 


20051) was f urther elaborated for the 


radia tion field of the Chakrabarti-Titarchuk type hybrid disc fjChakrabarti fc Titarchuk 


19951 ). Generally most of the papers which investigate the interaction of jet with disc pho¬ 


tons do not consider the issue of the launch mechanism of jets in terms of the accretion 


properties. In the advective disc regime, numerical simulations first showed t 


rat the extra- 


thermal gradient force in th e PSD automatically generates bipolar outflows flMolteni et al 


1994, 

1996; 

Das et. al. 

2014) 


201411 . Theoretically too, for a viscous advective disc, the mass out- 


flow rate of bipolar outflow was computed in terms o 


flChattopadhvav fc Das 


20071: 


Kumar 


fc Chattopadhva'^ 


various accretion disc parameters 


2013 


Kumar et al 


201811 . However, 


jets emerging from the PSD in the steady state and hydrodynamic limit are weak. On the 
other hand, if the emanati ng jets are simulta neously acted on by disc phot ons, then the 
l ets ob tained are stronger flKumar et r7,/.ll2014l ). In fact, it has been shown bv iKumar et al. 


(120141 ). that as the advective disc spectral state moves from low-hard or LH state to inter¬ 


mediate states or 


(Gallo et. al. 


2003 


M, the steady jet becomes stronger, as ha s been reported i n observations 


Fender et al. 


2010 


Rushton et al. 


2010l ). In other words. 


Kumar et al. 


(1201411 not only generated jets from accretion solutions, but also accelerated the jets by 
depositing m o mentu m of disc photons on to the jet. However, the formalism followed by 


Kumar et al. 


(120141 1 is only correct up to the first order of v/c, where v is the flow ve¬ 


locity. So solving the jet equations in the relativistic hydrodynamic limit is warranted for 
Ghakrabarti-Titarchuk type disc. 


(Hsieh & Spieerel 

1976; 

Mihalas & Mihalas 

19§4; 

Kato et al. 

1996) 


199811 . It was observed that in 


the non-relativistic limit, only flux of the radiation pushes matter. In the relativistic limit 
and optically thin plasma, the moving plasma is pushed by the flux and is dragged by the en¬ 
ergy density and pressure of the radiation field. The radiation transfers momentum on to the 
electrons and therefore, the radiative acceleration term is proportional to the number den¬ 
sity of electrons/positrons in the jet, which is true in relativistic, as well as, non-relativistic 
limit. But in relativistic limit, the radiative acceleration term is also proportional to the 
inverse of the enthalpy of the flow, which m arks a major d i fferen ce between relativistic and 


non-relativistic domain. Therefore, although 


Kumar et al. 


(120141 ) showed radiatively driven 


jets qualitatively explain the correlation between jet states and the spectral state of the ac- 
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Radiatively and thermally driven jets 5 

cretion discs, but a correct relativistic narrative is required. If however, one considers the gas 
pressure (p) to be negligible compared to the radiation pressure then even in the fully 

relativistic limit, the net radiative term is only proportional to the number density of the 
flow. Therefore, the physics of interaction of disc radiation with a fluid jet in the relativis¬ 
tic regime is qualitatively and quantitatively different than that be t ween radiation and jet 


in the non-relativistic domain (IChattopadhvav fc Chakra barti 


2002 - 


<umar et al 


relati vistic jets in the domain where Pij 3> p5ij (IChattopadhvav et al 


2004; 


2014), or 


ChattoDadhvav 


20051) . In this paper, we would like to investigate this phenomenon in details. 


Ferrari et al. 


(119851) studied interaction of fluid jets with the disc radiation but the jet fluid was described 
by hxed adiabatic index F and Newtonian gravitational potential. However, hxed F is an 
artifact of non-relativistic kinetic theory, i.e., if the internal ran dom motions of the particles 


in th e fluid is relativistic, then F is a function of temperature (jChandrasekhar 


1938 


Synge 


19571) . where, for relativistic temp eratur e s, the flow is described by F ~ 4/3 and for non- 

(119481 ) showed that it is unphysical to consider a 


relativistic temperatures F = 5/3. 


Tanb 


hxed F to describe huid with several orders of magnitude variation in temperature. Moreover, 


it has also been shown that, not only tempera 


depends also on the composition of t 


2009; 

Chattopadhvav & Chakrabarti 

2011 

Kumar & Chattopadhvav 

I 2 OI 4 

). Infact it 


ure, relativistic natu r e of the thermal energy 


le how (IChattoDadhvav 


Kumar et al 


2008 


2013 


Chattopadhvav &: Rvu 


Chattopadhvav et al 


201 


1 


Infact it was shown that, contrary to expectation, pair- 


plasma is thermally the least relativistic, and to make a how more relativistic one needs 
baryons in addition to electrons or positrons, and the huid with proton number density 
around 27% of that of electrons is thermally the most relativistic. Now radiation force is 
imparted mainly onto the electrons and positro ns which should make the lighter jet to move 


faste r due to inertia, as is shown for cold jets (IChattopadhvav et al. 


20041: 


Chattopadhvav 


20051 ). On the contrary, the thermodynamics of the jet will make the huid with 27% pro¬ 


ton content to be more relativistic. So in presence of radiation driving, would a jet with 
protons less than 27% be more accelerated, or, the thermal nature of the how dictate the 
dynamics. Additionally for huid jets, the pressure gradient force should accelerate the jet, 
while gravity decelerate, and radiative force depends on the huid speed in addition to the 
radiative moments. At the base, the jet is supposed to be hot, so initially, can the jet be 
thermally driven to speeds above Ueq, and therefore, is mostly decelerated by radiation and 
not accelerated at all? In short, can disc radiation power jets to relativistic terminal speeds? 
These are few of the various issues we would like to investigate in this paper. In this paper 
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6 Vyas et al. 

the accretion disc plays an auxiliary role, i.e., we do not employ a self consistent calculation 
of accretion-ejection system, rather the jet is assumed to originate in the inner disc, which 
in our case is PSD, and the accretion disc just provides the radiation held which interacts 
with the jet. While doing so, we have not increased the number of accretion disc parameters, 
and the shock location is estimated from the SKD accretion rate, and PSD luminosity has 
been estimated from the accretion rates from SKD and KD. 

In the next section, we present the simplifying assumptions and various equations. In 
sectionwe present the governing equations of the jet. In sectionwe present the model 
of accretion disc, estimation of the disc intensity from various disc components (section 
12.2. ip . and computation of the radiative moments fsection l2.2.2p from the disc intensities. 
Disc intensities, the disc velocity and temperature prohles are estimated in Appendix 
The post-shock intensity is estimated in Appendix [Bl In section |3l we present the solution 
methodology and in section 01 we present the results and in the last section ([5]) we discuss 
the physical implications of the obtained results and draw concluding remarks. 

2 ASSUMPTIONS, GOVERNING EQUATIONS OF JET & STRUCTURE 
OF ACCRETION DISC 

The space time metric is given by 

ds^ = —c^dt^ + dr'^ + r'^dcj)^ + dz^ ( 1 ) 

where, f, r, 0, 2 ; are the time, radial, azimuthal and the axial coordinates. The jet is 
considered to be in steady state {i.e., d/dt = 0). The above is a special relativistic metric. 
Although the jets originate close to the central compact object, but it traverses to distances 
where the effect of gravity is negligible. Hence pseudo potential is used to take care of gravity 
in the equations of motion, such that at close ranges the gravity limits the outward thrust, 
and at large distances it is just special relativistic regime, therefore one may avoid general 
relativistic complications and still be accurate enough. Moreover, the jets are collimated, 
so we consider on-axis {i.e., = d/dr = 0) and axis-symmetric {d/dcf) = 0) jet. 

Without any loss of generality, the jet is assumed to expand radially along the z-axis, the 
value of any jet variable on any particular point of the axis, is assumed to be maintained 
along its breadth at the same In the present effort, we just assume that there are jets 
from an accretion disc, but do not compute jets self consistently from accretion solution. 
The accretion disc is present but plays a supportive role, by supplying the radiation which 
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[h] 


SKD 



PSD 



\ Bipolar 
spherical 
outflow 

Figure 1. Cartoon diagram of cross-sections of axis-symmetric accretion disc and the associated jet in (r, (p,z coordinates). 
The Keplerian disc (KD), sub-Keplerian disc (SKD), and the post-shock disc (PSD) are broadly the three regions of the disc. 
The shock location Xs^ the intercept of SKD on the axis (do), height of the shock Hs are shown, the outer edge of the disc 
ro are all marked. Semi-vertical angle of PSD is 9ps and for SKD it is 0sk- The gradient of colour represents low temperature 
(red) to high temperature (blue). 


drives the jet. In the next subsection we present the equations of motion. We present the 
description of the accretion disc and the method to compute radiative moments at section 

o 


2.1 Equations of motion of the jet 

The energy momentum tensor of the jet and the radiation held is given by 

= (e + Tf = JRrfdudn, (2) 

where, suffix M stands for jet material and R stands for radiation held. The internal energy 
of the jet is e and p is the isotropic pressure of the jet huid, the metric tensor components 
are given by and m" is the four velocity. Furthermore, R is the specihc intensity of the 
radiation held, l°‘s are the directional derivatives, u the frequency of the radiation and is 
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8 Vyas et al. 

the solid angle subtended by the held point on to the source point. By dehnition, held point 
is where the moments are computed, while source point is the location of the source of the 
radiation. The equations of motion are given by 


Tf = 0 = + TtfU ( 3 ) 

From the above equation, the momentum conservation equation, in the direction is ob¬ 
tained by using projection tensor (gl^ -|- u^Ua), i.e., 

(g‘„ + u‘u,)T;1 = -(si + (4) 

Similarly, the energy conservation equation is obtained by taking 


rpa0 _ rjlO 

^a-L M.p — R-p 


-10.(3 


( 5 ) 


The derivation of the equations of motion of radiation hydrodynamics for optically thin 
plasma, using above preliminaries, was investigated by a number of workers. Since we study 
on axis jet, the equations of motion greatly simplihes. The momentum balance equation 
(equation 0 ]) , in steady state and for on axis jet becomes; 


{e + p) \u — -h 


GM, 


B 


dp z zdp , 

dz dz rripC 


dz ' [z — TgY^ 

The term containing Mg in the r. h. s of equation 




( 6 ) 


the gravity of non-rotating BH (jPaczvhski fc Wiita 
( 1 ^ in the scattering regime is, 
de e + pdp ^ 


is th e Paczyhski-Wiita term mimicking 


198011 . The energy conservation equation 


( 7 ) 


dz p dz 

where, p is the total mass density, is the leptonic mass density of the flow and irie is the 
electron rest mass. Similarly, from continuity equation the mass outflow rate is given as 


Mout = pR^A, A (X z^ for radial, narrow jet about the axis! 


( 8 ) 


In above equations, G, Mb, ct, Tg = 2 GMb/A‘ and A are the universal gravitational 
constant, the mass of the central black hole, Thomson scattering cross-section, Schwarzschild 
radius and cross section of the jet respectively. is the net radiative contribution and is 
given by; 


rr 

(Jj' ^ 


m c 


CTj' 

m. 


y- J U, J-l UjgJ. |U,1^ U,jJ 

c \ c 


- - UjP‘> + m' ( 


( 9 ) 


= [ 7(1 + 2Muz)J^ — j^R^S — V{uz + MuzRz)] for on axis jet! 
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Radiatively and thermally driven jets 9 

In above equations, E, and are the radiative energy density, the radiative flux 
and the radiative pressure tensor measured in observer frame, and S = —E, E = -^E^, 
and V = Furthermore, 7 (= —Ut = \/l + UiE) is the Lorentz factor. 

It may be noted that, we have assumed the jet to be flowing radially out within a conical 
surface for simplicity, since our primary concern is to investigate the respective role played by 
thermal and radiative driving terms on jet dynamics. Needless to say depending upon initial 
condition and disc radiation held, the jet geometry may depart from the simple geometry 
we are following in this paper. Even then, our assumption is not completely outlandish. The 
funnel like surface of the PSD is the region from where the jet is supposed to originate, the 
shape itself will arrest the lateral spread of the outflowing matter. Moreover, such shape 
causes the r component of radiative flux directed towards the axis (jChattopadhvavl 1200511 
which would also reduce the lateral spreading even with high jet-base temperature. However, 
one can justify our assumption only if we solve the jet equations since we need some estimate 
of the pressure. In appendix (JD]), we estimated the pressure gradient term along with the 
radiative term along r and compared them, the assumption of conical jet cross-section seems 
to hold. 


2.1.1 Equation of state and the final form of equations of motion 

The physics of the jet propagating in the radiation held of the accretion disc can be un¬ 
derstood, if equations (EHH]) are simultaneously solved. However, one has also to supply a 

bS) in order to 


closure relation i.e., a relation between e, p, p called the equation of state 


solve equations (EMI). An EoS 
(2008); Chattopadhvav fc Rvu 


or mul tispecies, relativistic flow proposed by 
(1200911 is adopted, and is given by. 


Chattopadhvav 


e = n^-meC^f 

with Ug- is the electron density and / is given by 


/ = ( 2-0 


1 + 0 


/90 + 3 

V30 + 2;j 




+ 0 


'90 + 3/7' 


( 10 ) 


( 11 ) 


7 \30 + 2/7^ 

Here, non-dimensional temperature is defined as 0 = fcT/(meC^), k is the Boltzmann con¬ 
stant and f = Hp+fue- is the relative proportion of protons with respect to the number 
density of electrons. The mass ratio of electron and proton is 7 = mpjmp +. It is easy to se e 
that by putting .^ = 0, we generate EoS for relativistic e~ — e+ plasma flRvu et a/.l 1200611 . 
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10 Vyas et al. 

The expressions of the polytropic index N, adiabatic index T and adiabatic sonnd speed a 
are given by 


N = 


Idf 


r = 1 + 


Tp 


2r0 


( 12 ) 


2de’ ' ' ' N’ e + p f + 2e 

This EoS is an approximated one, and the comparison with the exact one shows that this 
EoS is very accnrate (appendix [C]). Additionally, being algebraic and avoiding the pres¬ 
ence of complicated special fnnctions, this EoS is very easy to be implemented in simn- 


lation codes, as well as, be nsed in analytical investigations lIChattoDadhvav fc Rvn 


Chattopadhvav fc Chakrabarti 


2011 


Rvn et al. 


2006 


Chattopadhvav et al. 


2009; 


20131) . The jet 


plasma is fnlly ionized. Therefore the interaction with photons wonld be dominated by scat¬ 
tering. Therefore, the energy eqnation ([7]) has no sonrce term becanse in the scattering 
regime and in absence of emission/absorption, the r. h. s is zero and the flow is isentropic 
flMihalas fc MihalasI ll984J ). Under snch conditions, eqnation ([7]) along with eqnation flTOl) 
can be integrated to obtain the relativistic isentropic eqnation of state, 

p = Cexp(fc3)03/2(30 + 2)^1 (30 + 

where, ki = 3(2 —^)/4, ^2 = 3^/4, = (/ —r)/(20), r = (2 — ^ + ^/h) ^-nd C is the constant 

of entropy. We replace p from the above eqnation on to eqnation ([8]), we get the expression 
for entropy-accretion rate, 

Mont 




= exp(fc3)03/2('30 ^ 2)^1 (30 + 2/p)^^u^z^ 


(13) 


geom. const. C 

This is also a measnre of entropy of the jet and remains constant along the jet. We adopt 
a nnit system where, the nnit of speed is c, nnit of length Vg = 2 GMb/c^ and the nnit 
of mass is Mb- Henceforth we write all eqnations in this nnit system, except where it is 
explicitly mentioned. The three-velocity v is given by = —UiU^/utv} = —UzU^/utu'', i.e., 
= Uz = 'jv. Now nsing energy conservation eqnation ((71) along with the eqnation of state 
CB. the gradient of temperatnre of jet is given by. 


dQ 

dz 


0 

'n 


v\dz z 


The momentnm balance eqnation 
dv 2'y^c? 


(14) 


), with the help of eqnations (|T0|1 , (|T2|l and flTTl) , becomes 


4 

7 V 


\--A- = 




dz 


+ 


^ 2(^-l)2 

7^(2-Or., , . 


/ + 20 

= a^: -|- a^ -f- a^. 


[(l + v^)F -v[£ + V)] 


(15) 
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Radiatively and thermally driven jets 11 

The 1. h. s is the net acceleration term of a steady state jet. On the r. h. s, the hrst 
term is the thermal term a^ = ‘I'y^a^ jz and it accelerates, while the second being gravity 
Big = —0.5/(z — 1)^, it decelerates. The third term in r.h. s is the radiative term a^ = 
7^r[(l + v‘^)J^ — v(S + V )]/(/ + 20). The radiative contribution is within the square bracket 
and the rest represents the interaction of matter jet with the radiation field. The physical 
signihcance of the term in the square bracket term is worth noticing. It has the form 

{l + v^)J^-v{£ + V) 


The term proportional to v comes with a negative sign and would decelerate and is called the 
radiation drag term. If the first term (1 + v‘^)jF dominates, then radiation would accelerate 
the flow, which means the net radiative term would either be accelerating or decelerating 
depending on the velocity. The dependence of radiative term on v arises purely due to 
relativity. In the purely non-relativistic domain i.e., n 1, the radiative term is just jF. In 
the fast but sub relativisti c domain i.e., v'^ 1 the radiative term is jF — v(S + V) sirn ilar to 


the formalism followed by IChattopadhvav fc Chakrabartil (1200211 ; 


Kumar et al. 


(120141 . The 


drag term arises due to the resistance faced by the moving material through the radiation 
field, and the finite value of the speed of light. Much talked about equilibrium speed Ueq is 
when a.r = 0, i.e., 

S + V 


-"eq 


= 3? — \/3? — 1; where, = 


2J^ 


(16) 


From equation (ITOi) . it is clear that if the relative contribution of radiative moments or 3? 
approaches 1, i.e., T = E = V, then Ueq —)■ 1, i.e., no radiation drag. Therefore, the nature 
of the quantity 3? dictates, whether a radiation field would accelerate a flow or decelerate it. 
Of course the resultant acceleration depends on the magnitude of all moments. There is an 
added feature of radiatively driven relativistic fluid, i.e., the radiative term is multiplied by 
a term inverse of enthalpy ({/ + 20}/r) of the flow, which actually suggests that the effect 
of radiation on the jet is less for hotter flow. 


2.2 Accretion disc and radiative moments 


The acc retion disc model considered here is the hybrid disc of th e Chakrabarti-Titarchuk 


flavour (jChakrabarti fc Titarchuk 


1995; 


Giri fc Chakrabarti 


291311 . In Fig. [H we show all 


the components of the disc, i.e., PSD, SKD and KD. The SKD flanks the KD, but mingles 
and forms the single component PSD. The colour coding represents lower (red) to higher 
(blue) temperature. The outer edge of the disc is tq where the disc height is Hq. The inner 
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12 Vyas et al. 

edge of SKD and KD is the shock location Xg, the inner edge of PSD say is in principle 
the horizon, bnt we have considered it to be ri^ = while calcnlating the radiative 


moments, since very little radiation is expec 


shoe 


i height is marked as Hg = 0.6 (x^ — 1) fjChakrabarti fc Titarchnk 


ed from a region very close to t 


1995 


le horizon. The 


Chattopadhvav 


20051) . therefore ^ps = tan ^{xg/Hg). The semi vertical angle of the SKD (^sk in Fig[T]) is 


taken to be 85°. This assnmption has been dictated by a lar ge nnmber of simnl a tions 

1 


which showed SKD to have a flatter surface compared to PSD (iMolteni et al. 


199 ^ 


Giri & Chakrabarti 

2013; 

Das et. al. 

2014) 


19961: 


20141) . The intercept of the SKD surface on the z- 


axis do = 0.4 x Hg. We chose ro=3500rg. SKD emits via synchrotron and bremsstrahlung 
processes, so the information of velocity {i.e., density) and temperature prohle is required, 
and have been estimated in appendix A. Injection speed {'do) for the SKD at the outer disc 
boundary is kept 0.001, the angular momentum of the disc A is 1.7 and the temperature at 
r = To is ©0 = 0.1. The PSD is hotter than the rest of the dis c (including SKD and KD ) and 
puffs up in the form of a torus. KD emits thermal photons fjShaknra fc SnnyaevI 119731 ) and 
SKD emits via bremsstrahlung and synchrotron processes. PSD emits bremsstrahlung and 
synchrotron photons, as well as, being fatter and hotter, inverse Comptonize these photons 
and the photons intercepted from SKD and KD to produce hard radiation. All the spectral 
states therefore, can be obtained by controlling the SKD accretion rate Mgk and the KD 
accretion rate Mkd. If Mkd is relatively less than Mgk, then due to the lack of supply of 
soft photons, PSD will remain hot, and thus producing the low hard or LH state. Increase 
in viscosity and/or increase in accretion rate at the outer bounda ry would push the shock 


closer to the central object (jKumar fc Chattopadhvavll2oi^. l2014f) . This would brighten up 


the disc, but would make the spectra softer as the number of hard photons from the PSD 
would be lower due to the decrease in PSD size. Hence the spectral index would increase 
and spectra would enter the intermediate or IM states. The increased size of the pre-shock 
disc (SKD-I-KD) and the increased accretion rate, would eventually weaken the PSD or 
completely destroy it, the contribution of hard photons would plummet and the disc would 
be more luminous, while the spectra would become soft similar to a multicoloured black 
body. 
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Radiatively and thermally driven jets 13 


2.2.1 Relativistic transformations of intensities from various disc components 


In order to compute radiative moments, we need to k now radiative intensities o f various disc 


components. The intrinsic KD intensity is given by flShakura fc Snnvaev 


19781) 


.^kdO — 





erg cm ^s ^ 


(17) 


To compute the radiative moments from SKD, we need to know the temperature and 
density distribution of SKD, in order to calculate the intrinsic radiative intensity of SKD. 
The density and temperature of SKD starting with some outer boundary condition can only 
be solved numerically. However for simplicity, we estimate the approximate values of velocity, 
density and temperature prohle of SKD in order to compute the Jsko- Equations flAll IA3jl 
give the analytical expression of all the components of three-velocities and the corresponding 
Lorentz factors (Appendix of the SKD. This allows us to compute the density prohle for 
a given Mgk. The density and temperature prohle of SKD are estimated in equations flAbI fc 
IA7I1 in Appendix lAll 

We further assume that there is stochastic magnetic held in the SKD which is in partial 
equipartition with the gas pressure. The ratio of magnetic pressure (pmag) and the gas pres¬ 
sure (pgas) is also assumed to be constant fd i.e., Pmag = = (dpgas = /Snsk^Tsk, where, 

nsk and Tgk are the SKD local number density and temperature, respectively. The emission 
mechanism is dominated b y synchrotron and bremsstrahlung emission, and therefore the 


SKD intensitv is given bv 

f Svensson 

1982; 

Shaniro & Tenkolskv 

1983; 

Kumar et al. 

2014 

Kumar & Chattonadhvav 

2014 

) 


-^sko R-vn T h 


syn 


‘brem 


16 

/eHsk\ 

3 c 

V rUeC J 


^- 27,^2 


'©skrUe 


(do sinOs-k + r cosOsk) 


erg cm ^s ^ 


(18) 


where, ©sk^nskjT, 6*sk, do, B^k and gi, are the pre-shock local dimensionless temperature, 
electron number density, horizontal distance from center of the disc, angle from the axis of 
symmetry to the pre-shock surface, the intercept of the SKD surface on to the axis of sym¬ 
metry, the magnetic held and relativistic Gaunt factor {gh = 1-|-1.780^'^^), respectively. The 
factor outside square brackets serves as the conversion factor from emissivity (erg.cm“^s“^) 
into intensity (erg.cm“^s“^). 

The expression of intensity is more complicated for PSD, as a result we make further 
simplifying assumptions. The PSD itself emits via bremsstrahlung and synchrotron, but also 
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14 Vyas et al. 

inverse Comptonizes its own photons, as well as, photons intercepted from SKD and KD. It 
is beyond the scope of this pape r to do a proper radiative trans fer treatment of the accretion 


disc. Instead, w e fed the code of 


five solutions of 


Mandal fc Chakrabartil (120081) with the viscous and dissipa- 


Kumar fc Chattopadhvavl (1201411 as back ground solution, and computed the 


spectra from PSD, SKD and KD. The shock location Xs is estimated from equation IBII and 
is also presented in Fig. |B])a). Although the relation between Xg a nd Msk has been obtained 
by ge neralizing the solutions for a certain viscosity parameter (iKumar fc Chattopadhvav 
20141) . but we treat them as generic, since the general pattern is similar for a large number 
of cases we have analyzed. Then we fit the ratio of the luminosities of PSD and that from 
the pre-shock disc (SKD and KD) y, as a function of Xg from the spectra obtained by the 
radiative transfer of the background solutions. The relation between y and Xg is given by 
equations IB41 IB 3 1 and is also plotted in Fig. [B](b). The preshock luminosities are obtained 
by integrating 4^ and Ad over the respective disc surfaces. So in principle we have two free 
parameters to fix the radiation held above the accretion disc, ilAk and Mkd. The intensity 
as measured in local rest frame of PSD is given by 


f^ps/^Aps As.f^Edd/^Aps (ergcm s ), 


(19) 


where Lpg and Apg are the PSD luminosity and the surface area of the PSD respectively. 
L^dd is the Eddington luminosity and £ps is the PSD luminosity in units of pEdd- 

The intensities from either the PSD or SKD or KD components of the disc, namely, 
Ipso, Ako and AdOi are measured in the local rest frames of the disc components. However, 
the matter in the disc is moving, so one has to transform the quantities into the observer 
frame! The intensity measured in the observer frame is presented in compact notation as. 


= A 


‘■JO 


( 20 ) 


V [1 + 

Here 7j is Lorentz factor and v® is component of 3-velocity of accreting matter and /®s 
are directional cosines. The suffix j —)■ ps, sk, kd signihes the contribution from either PSD, 
or SKD, or KD. For PSD and SKD v® is calculated following Appendix lAj while for KD 
V® = (0, Vkd, 0) is the Keplerian azimuthal velocity. The luminosity from various components 
of the disc is obtained by integrating the respective local specihc intensities over the disc 
surface i.e., 


L ,=2 


Ij X 27rr cosec ^6*j 


dr. 


( 21 ) 
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Radiatively and thermally driven jets 15 

here, j represents various disc components. To compute the luminosity, the limits of inte¬ 
gration are inner and outer limits of the disc components, for example for PSD the limit of 
integration is rjn —>■ while for SKD and KD the integration limits are Xg —)■ Tq. Apart 
from integration limits, various disc components are identified by the respective 0jS, which 
defines the surface of disc components. The total luminosity is given by L = Lpg-I-Tsk + Tkd, 
and in units of Eddington limit it is £ = £ps -1- 4k + 4d- All the transformations presented 
above, are exactly correct in the special relativistic regime. We are not taking into consid¬ 
eration the phenomenon of light bending since the jet spans from close to the horizon to 
few thousands of Vg. Beyond few tens of Vg the general relativistic effect may be ignored, 
but special relativistic effects cannot be. But close to the horizon, if only special relativ istic 
effects are considered then the Jps gets unnecessarily jacked up fIChattopadhvavl l2005l) . In 
order to address this, we estimate the fraction of radiation that will not be absorbed by the 
black hole and escape. From geodesic equations of p hotons, it can be easily sho wn that if 
smip > 3\/3(l — \ jw') jijlw') then the photon escapes (IShapiro fc Tenkolskvlll988l ). where 4 
is the angle of the direction of propagation of light with the radial direction, and the w is 
the spherical radius coordinate. So we express 4 in terms of r and ^ps the semi vertical angle 
of the PSD inner surface. Assuming locally the radiation is isotropic, then the fraction of 
intensity from PSD which would escape and interact with the jet is, 

TT — sin ^ (^3\/3(sm6'pg/2r)(l — sinO^g/r)^ 


n = 


TT 


( 22 ) 


Since SKD and KD is further away from the black hole, no such estimation of gravity effect 
on emitted radiation from these components is needed. 


2.2.2 Computation of radiative moments 

Radiative moments are the zeroth. Erst and second moments of specihc intensity, and fre¬ 
quency integrated moments are respectively called radiation energy density, radiative flux 


and radiation pressure, and are expressed as following. 


^ E ^ 


( i \ 

Q 


Fi 

-II 

r 


V ^ ) 


\ — / 


lydvdVL 


Here, /*s are directional derivatives, v the frequency of radiation and D is the solid angle 
subtended by the field point (where the moments are calculated) on to the source of radiation. 
Since we are considering on axis jet, therefore we need to compute the radiative moments 
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16 Vyas et al. 

only along the axis. The radiative moments along the jet axis {z axis) are calculated from the 
PSD, SKD and KD components of this hybrid disc model, these components are indicated 
in following expressions with subscripts ps, sk and kd respectively. 


^ ~ nw (y J + J Ikddi^kd^ — — {Eps + T^sk + .^kd) 

= + ^sk + i’kd (23) 

^ = — ( / + [ Isk^dDsk + / /kd/"df2kd) = — (Fps + Fsk + Fkd) 

= -Tps + J^sk + -Tkd (24) 

v= — (f Jps/^/^dDps + / IskEEdQ,^ + f hJ^Edn^^) = ^ (Pp^ + Psk + Pkd) 

TTIjC J J / TTIj 

= Eps + Esk + Pkd (25) 


All the points on the axis of symmetry are held points {z), i.e., where radiative moments 
are to be computed. The coordinates on the disc surface are r, 0, z' and z' = r cot6*j, where 
j = ps/sk/kd. It is easy to see that for extended source, and held point close to the source, 
the directional cosines of E < 1, but for point source = 1, and E = T = V. Therefore, 
according to the note following equation (1T5|1 . for point sources 3? = (£^ + P)/2P = 1 and 
there will be no radiation drag. The radiative moments from PSD are, 

TZzrdrdcj) 


r2-K 


Eps — S 


Pps = S 


0 [{z-r COtOpsY + r2]3/2y4^ [1 + Vj/*] 

2^^ 7lz{z — r cot6ps)rd(pdr 


4 

ps 


I 


Pps — *5 


0 [{z-r COtOpsY + r2]2y4^ [1 + Vi/*]p^ 

xs /•27r TZz{z — r cot6psYrd(j)dr 


! i 

Jri„ Jo 


( 26 ) 

(27) 

(28) 


[(j - r coMp.)2 + (1 + v,;‘)J, 

iS is a constant, which is obtained after converting the whole expression in geometric units 
mentioned after equation | 


^ 1.3xl03%,a^ 

27rcmeApsGMQ 

ax, me, G, Mq, Apg, ips are Thomson scattering cross section, rest mass of electron, constant 
of gravitation, solar mass, surface area of PSD and post shock luminosity in units of Edding¬ 
ton luminosity (pEdd = 1.3 x 10^^Mb/Mq) respectively. Now, depending on the central mass, 
ips is calculated from equation flB3l or [B4l) . and hence depends on SKD and KD accretion 


rates. 

The radiative moments from the SKD are 
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— 


= 


Vak — 


With 


-^sk ‘^ks 


^ {rcosOsk + dosinOskjzdrdj) 

lo {cote^kr + d^f [{z - r cote^^^f + r^f/^'jX (1 + v*/*)! 

/■2’r ^ {rcosOsk + dosin9sk)z{z — r cot9sk)drd(j) 

Jo {cot9skr + dof [{z - r cot9sk)‘^ + (1 + 

^ {rcos9sk + dosin9sk)z{z — r cot9sk)‘^drd(j) 

Jo {cot9^],r + d^f [{z - r cot9^],Y + (1 + 


fUoroHo\'^^^^ [l M (UoroHo^^"^ 


Constants iSks and 5kb are associated with synchrotron and bremsstrahlnng terms and are 


given as 


iSks — 


9.22 X lO33e403/3arm,\ 
TTm'lm'^+c^G'^ Mq 


iSkh = 


1.51 X lO^mg]^\/0^(TT 

n'^mp+MQG'^y/mek 


Here rhsk = Mgk/MEdd, where the Eddington accretion rate MEdd = 1-44 x 10 ^’^(Mb/Mq) g 
s“^ and hEdd = ^EddC^- From a certain point x on the axis, due to the shadow effect of the 

PSD the inner edge of SKD observed is given by 

. X __ z -dp _ 

^ {z - Hs)/Xs + COt9sE 

Therefore, to obtain the radiative moments at some x on the jet axis the integration limit 


on r is to rp. 

Similarly from KD, the moments are 
ro /-aTT z{r~‘^ — \/3r~^/‘^)d(pdr 


Tkd — /C 


JTd = 1C 


Pkd =/c 


Jo {z^ + (1 + 

n 27r z{r~'^ — y/3r~^^‘^)d(j)dr 
0 (^2 + r2)274^ (1 + 

n2-K z(r~^ — \/3r~^^‘^)d(j)dr 

0 (x2+r2)5/2^4^(l + v,/*)J 


The shadow effect of blocking a fraction of radiation from the KD by the PSD is also 
taken into account and the integration is done from r\^ and is given by 


and the dimensionless constant /C is given by 

4.32x lO^^rhkdO'TC 

327r2meGMQ ’ 

with rhkd being the Keplerian accretion rate in units of MEdd- 
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Figure 2. Distribution of radiative moments from the post-shock disc (PSD) with 2 : (in units of Vg). Each curve represents 
^ps/S (solid, black online), J^ps/S (dotted, red online) and Vps/S (dashed, blue online), respectively. Various panels are for 
radiative moments (a) without relativistic transformations in the accretion disc, (b) with special relativistic transformations 
up to first order in u, (c) with full special relativistic transformations and (d) with approximate general relativistic correction 
and full special relativistic transformations. All the figures are obtained for rhgk = 5.06, which produces a shock at Xg = 20 
around a IOMq BH. 


2.2.3 Nature of radiative moments 


We now numerically integrate equations (12611351) to obtain the radiative moments from PSD, 
SKD, and KD. In Fig. |2li-d we present a comparative study of the space dependent part of 
the radiative moments i.e., EpsIS (solid, black online), J^ps/S (dotted, red online) and VpsjS 
(dashed, blue online) from PSD which were computed in previous works and the present one, 
because we would like to show the effect of various corrections considered, while computing 
the radiative moments. To compare the moments from PSD, we need to know the size of PSD. 
This is obtained by estimating the Xg from the supplied rhsk- For rhsk = 5.06, we estimate 
Xg = 20 (equation IBip . Figure [2]( a) presents the radiative moments from PSD when the 


accre tion disc was treated as a simple radiator without any dynamics (jChattopadhvay et al. 


2004J ). In Fig. Efb), we applied special relativistic transformations while computing the 


intensities in the observer frame from the local disc frame, but the Lorentz factor (7^) in 
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Figure 3. Distribution of radiative moments S (solid, online black), T (dotted, online red) and V (dashed, online blue) from 
PSD for (a) IOMq and (b) 10® Mq black holes. Distribution of radiative moments from the Kd (c) and from SKD (d) is same 
for both types of black holes when expressed in the geometric units. Various parameters used to compute the moments are 
rhgk = 7, rhkd = 1 and (3 = 0.5. This produce Xs = 13.2032 and luminosities are = 0.0265, £kd = 0.039 and for stellar mass 
BH ^ps = 0.215 (a), while for larger BH ^pg = 0.661 (b). 


the red-shift factor of the intensity (equation [20]) was ignored fIChattopadhvavl l2005l) . And 
in Fig. [2ic) we show moments computed with correct special relativistic transformations. 
Figure [2id) shows only the distribution of space dependent part of radiative moments from 
PSD which is being used in the current study. In this panel, moments are with full special 
relativistic transformations and general relativistic correction for the radiation absorbed by 
the black hole near horizon. We notice that general relativistic correction (Fig. [2li) reduces 
the absolute value of the moments and the moments peak at a slightly different location on 
the axis than the previous case (Figj2j::). Ignoring the effect of disc motion while computing 
radiative moments, under-estimates the moments (Fig. [2^), while partial implementation of 
the very effect, over estimates the radiative moments (Fig. 12b) to unrealistic values. In this 
paper, we use full special relativistic transformation and general relativistic corrections for 
the radiative moments from PSD ( as in Fig. [2li). 
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Figure 4. Combined moments from (a) 10 and (b) IO^Mq black holes. Various curves are £ (solid, online black), (dashed, 
online red) and V (dashed, online blue). Disc parameters are mgk = 7 and rhkd = = 0.5 for which the shock obtained is at 

Xs = 13.2032, luminosities are = 0.0265, — 0.039. The £ps = 0.215 for IOMq BH and ips = 0.661 for 10 ^Mq BH. The 

moments are expressed in geometric units. 


From Appendix IB 11 it is clear that £ps is different for IOMq and IO^Mq BH, for the same 
set of free parameters i.e., rhsk and rhkd- This should affect the net radiation field above the 
disc. In Fig. |3l,, we plot Tps (solid, black online), J^ps (dotted, red online) and V^s (dashed, 
blue online) with z for rhsk = 7 and rhkd = 1 for Mb = IOMq. For such accretion rate the 
shock is at Xs = 13.203 (see, equation IBip . In Fig. Eb, we plot radiative moments above a 
disc around a Mb = 10®Mq BH, for the same set of accretion parameters. The radiative 
moments from the PSD around 10 ®Mq BH are about three times than those around IOMq. 
It may be noted, that for lower rhsk the shocks are formed at larger distance away from the 
BH, and the radiative moments around stellar mass and super-massive BH are same. In Fig. 
El, we plot £^kd (solid, black online), J^kd (dotted, red online) and Pkd (dashed, blue online). 
In Fig. Eb we present Tsk (solid, black‘ online), J^sk (dotted, red online) and Vsk (dashed, 
blue online). The PSD luminosity for IOMq BH is about f'ps = 0.215 and f'ps = 0.661 for 
IO^Mq BH. The luminosities of the pre-shock disc 4k = 0.0265 and 4d = 0.039 are same for 
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discs around super massive, as well as, stellar mass BH. The moments due to PSD around 
a super-massive BH is larger than that around stellar mass BH, however, the SKD and KD 
contributions in the geometric units are exactly same for stellar mass and super massive 
BH. In physical units these moments would scale with the central mass. Finally, we show 
combined radiative moments from all the disc components for IOMq (Fig. 0^) and IO^Mq 
BH (Fig. |3 )d) for exactly the same disc parameters as in Fig.[3K-d. For higher rhsk, the overall 
radiation field (in geometric units), above a disc around a stellar mass BH is different than 
the moments around a super massive BH. This is because the for higher rhsk the shock in 
accretion is located closer to the BH, which produces a cooler PSD around a stellar mass 
BH than a super massive BH and therefore, larger efficiency of Comptonization. However, 
for lower rhsk (equation [Bl] & Fig. lAlh ). the shock is located at larger distance from the 
BH, and the efficiency of Comptonization is similar for both kinds of BH, and hence the 
moments are similar too. 


3 SOLUTION METHOD 

To obtain the solution of steady state, relativistic jet, we need to integrate equations flTTl) 
and f[T5]l simultaneously. Alternatively, one may integrate equation flT^ with the help of 
equation o. since equation flT^ is the integrated version of equation ffTT|) . We employ 
Runge-Kutta’s 4*^ order method to integrate differential equations. 


3.1 Sonic point conditions 

Since jets originate from a region in the accretion disc, close to the central object, the 
base jet velocity should be small. However due to this proximity of the jet base to the 
central object, the temperatures at the jet base should be very high. In other words, jets 
are subsonic at its base. While far away from the central object the thermal energy and 
the radiative energy would drive jets to large v, but simultaneously becoming less hot, i.e., 
outer boundary condition of the jet is super sonic. This means at some point the jet would 
become transonic, and the point in which this happens is called the sonic point Zc and the 
derivative dv/dz —)■ 0/0 (eq. [15]). This gives us the so-called sonic point conditions. 


Vr — ttr 


a. = 


27 / 


7/^ 


2(Ze-l)2 /c + 20 


^ {(1 + ~ Vc{Sc + 'Pc)} 


(36) 
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Figure 5. Variation of ac (a, b), ©c (c, d) and A4c (e, f) as a function of Zc for jets are around Ml = IOMq BH (a, c, e) 
and M8 = IO^Mq BH (b, d, f). Each curve corresponds to rhgk = 13 (solid, red online), 10 (long dashed, magenta online), 8 
(dashed-dotted, green online), and only thermally driven jet (long-dashed dotted, black online). For all the plots = 1. 


The dv/dz\c is calculated by employing the L’Hospital’s rule at Zc and solving the result¬ 
ing quadratic equation of dv/dz\c- The resulting quadratic equation can admit two complex 
roots leading to the so-called O type or ‘centre’ type sonic points, or two real roots. The 
solutions with two real roots but with opposite signs are called X or ‘saddle’ type sonic 
points, while real roots with same sign produces the nodal type sonic point. The jet solu¬ 
tions flowing through X type sonic points are physical, and in this paper care has been taken 
to study jet solutions through X type sonic points. So for a given set of flow variables at 
the jet base, a unique solution will pass through the sonic point determined by the entropy 
A4 of the flow. For given values of inner boundary condition i.e., at the jet base Zb, Vb and 
ttb, we integrate equation flT^ and ffTT|) . while checking for the sonic point conditions (equa¬ 
tions [36]). We iterate till the sonic point is obtained, and once it is obtained we continue to 
integrate outwards starting from the sonic point. From equation (l36jl it is clear that for a 
thermally driven jet, sonic point exists from Zc = 2 —)■ oo. However, radiatively driven flow 
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may not posses sonic point at large distances away from the jet base, because the presence 
of strong radiation held may render Oc 0 at those distances. In Figs. Eti'-f, we compare the 
how quantities Oc (a, b), ©c (c, d), Me (e, f) as a function of Zc- The left panels shows the 
sonic point properties of jets around IOMq BH (a, c, e) and the right panels show sonic 
point properties of jets around 10^Mq BH (b, d, f). The KD accretion rate, or, rhkd = 1 
is kept invariant for all these plots, but various curves are for rhsk = 13 (solid, red online), 
10 (long dashed, magenta online), 8 (dashed-dotted, red online), and only thermally driven 
jet (long-dashed dotted, black online). It is interesting to note that, the region outside the 
central object available for sonic points shrinks, as the disc luminosity increases. For lumi¬ 
nous discs say rhsk > 10, sonic points can only form for Zc < 8. This implies that only very 
hot how has thermal energy density comparable to the radiation pressure, and therefore for 
any how with less thermal energy may be considered as collection of particles rather than 
a huid in such radiation held. Moreover, for rhsk = 13, multiple sonic points may form for 
some values of Ai. 


4 RESULTS 

A radiatively inefficient disc can only give rise to thermally driven jets and not radiatively 
driven jets, so we choose luminous disc. We discuss the jet properties for electron-proton 
jets i.e., ^ = 1, until specihed otherwise. We choose rhkd = 1 to generate the KD radiative 
moments until specihed otherwise. The accreting material is assumed to posses stochastic 
magnetic held with constant magnetic to gas pressure ratio (3 = 0.5, until specihed otherwise. 
In Fig. [6^, we plot the jet 3-veIocity v (solid, black online) and the sound speed a (dashed 
dotted, blue online) as a function of This jet is from a disc around a stellar mass BH. The 
sonic point Zc is at the crossing point of v and a. The terminal speed achieved for this case is 
vt ~ 0.45, and the total disc luminosity is £ = f'ps -1- 4k + 4d ~ 0.14 in units of LEdd- In Fig. 
(Ob we plot AI as a function of z and in Fig. [61:: we plot the variation of F or adiabatic index of 
the jet. The inner boundary condition is Zh = l.br^, Vh = 0.014 and ah = 0.51, and the SKD 
and the KD accretion rates are rhsk = 5 and rhkd = 1- Since the interaction between radiation 
and the jet material is assumed to be in the Thompson scattering regime, the source term 
of the hrst law of thermodynamics turns out to be zero (he., equation flTll . and therefore 
AA. which is a measure of entropy, remains constant through out the flow. The base of the 
jet is very hot, therefore F —)■ 4/3 at the base. However, as the jet expands to relativistic 
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Figure 6. (a) The variation of jet 3-velocity v (solid, black online) and sound speed a (dashed-dotted, blue online) with 2 . 
(b) Ai is plotted as a function of 2 and (c) the variation of F with 2. All the plots are generated for inner boundary condition 
2 :^ = 1.5, Uf, = 0.014 and = 0.51, and the SKD and the KD accretion rates are rhsk = 5 and mkd = 0 = 0-5 and shock 

obtained is at Xs = 20.27rg and luminosities are £sk = 0.0082, £ps = 0.113, £kd = 0.027. (d) Comparison of v for a thermally 
driven jet (dotted, red online) and radiatively plus thermally driven jet (dashed, blue online). Here accretion parameters are 
rhkd = ^10 = 0.5 and rhsk = 10 for the radiatively and thermally driven jet. For thermally driven flow no radiation interaction 
has been considered. The luminosities are ^gk = 0.108, £kd = 0.055 and £ps = 0.38. The composition of the jet is ^ = 1 or 
e~ — p'^ plasma, and is launched around a IOMq BH. 





velocities (at 2; —)■ large), the temperature falls such that T —)■ 5/3. In Fig. |6]i, we compare 
the V prohle of a thermally driven jet (dotted, red online) with a radiatively plus thermally 
driven jet (dashed, blue online) starting with the same base values. The radiatively driven 
fluid jet (blue dashed) is powered by radiation from a disc with parameters rhsk = 10 and 
rhkd = 1- From the base to hrst few r^, the v prohles of the two flows are almost identical, 
and the radiative driving is perceptible at 2; > 7.66. The terminal speed of the thermally 
driven flow is slightly less than 0.45 and for the radiatively driven flow it is vt ~ 0.65. The 
radiative driving of the jet is ineffective in regions close to the Z},, because the thermal driving 
accelerates the jet to n ~ Ueq- Therefore radiative driving is ineffective in those region, and 
results in a similar v prohle up to z ~ 7.66, but beyond it radiative driving generates a how 
with a 44% increase in vt- 

In Fig. [T^-b, we show a transonic jet from a disc around a stellar mass BH, in which. 
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Figure 7. (a) Variation of v (solid, online black) and a (dashed-dotted, online blue). The base speed V}) = 0.087 and sound 
speed a}) = 0.545 at jet base = l.bvg, critical point obtained is at Zc = 2.^rg\ (b) Variation of at (solid, blue online), ar 
(dotted, black online) and a^ (dashed, red online) with 2 . The SKD and the KD accretion rates are rhsk = 12 and rh-kd = 1? 

= 0.5 and shock obtained is at Xs = 5.87 and luminosities of various disc components around IOMq BH, are = 0.295, 
.^pg = 0.522 and ^kd = 0.0667. The composition of the jet is ^ = 1. The vertical dashed line (magenta online) shows the position 
of sonic point Zc- 

the jet has been accelerated, as well as, decelerated by radiation. In Fig. [Tii, we plot v 
(solid, black online) and a (dotted dashed, blue online) as a function of ^ for inner boundary 
condition Vb = 0.087, ah = 0.545 at Zb = 1.5, and rhsk = 12 and rhkd = 1- In this case 
the sonic point is obtained at Zc = 2.5. The total luminosity turns out to be £ = 0.884. 
The 3-velocity v increases beyond Zc and up to 2; ~ 4, and then decelerates in the region 
^ < z ^7 and thereafter again accelerates till it reaches terminal value at vt ~ 0.86. Let us 
analyze the various terms that influence v. In Fig. [7 |d, we plot the variation of gravitational 
acceleration term or a^ (dashed, red online), the radiative term a^ (dotted, black online), 
and acceleration due to thermal driving a^ (solid, blue online). In both the panels the dashed 
vertical line (magenta online) shows the location of Zc- From the 1. h. s of equation ([15]), it 
is clear that in the subsonic region v can increase (he., jet accelerate) with z only if r. h. s 
is negative. While in the supersonic region the jet accelerates if the r. h. s is positive. The 
gravity term or a.g is always negative, while a^ is always positive. In this particular solution 
a^ < 0 for 2; < 8.53. In the sub sonic region he., z < Zc, |ar| <C a* and a* < a^, therefore r. h. 
s of equation flT^ is negative and the jet is accelerated. At the sonic point a* = a^ -|- a^. For 
© 0000 RAS, MNRAS 000 , 000-000 



I I I I I I 11 [ ' r 









26 Vyas et al. 



Z 




Figure 8. (a) Variation of v with 2 :. Each curve represents /3 = 1.0 (dotted, black online), 0.5 (solid, black online) and 0.0 
(dashed, black online), (b) Variation of with j3. Other parameters same as (a), (c) Variation of v with z. Each curve 
represents rhgk = 1.0 (dashed, blue online), 7 (dotted, blue online) and 10 (solid, black online), (d) Variation of ux with rhgk. 
Other parameters same as (c). (e) Variation of v with 2 . Each curve represents ^ = 1.0 (solid, black online), 0.6 (dotted, magenta 
online) and 0.3 (dashed, magenta online), (f) Variation of ux with Other parameters same as (e). Accretion parameters are 
rhgk = 10 (in a, b, e, f), {3 = 0.5 (in c, d, e, f), ^ = 1.0 (in a, b, c, d) and mkd = !• Jet base values are 25 = 1.5, U 5 = 0.014 and 
= 0.51. The mass of central black hole is IOMq. 


z > Zc, &r decreases to its minimum value at z = 5.15. Gravity is less important at these 
distances and a^ ^ at, which makes the r. h. s negative. Therefore, the jet decelerates in the 
range 4.78 < z < 6.34. For z > 6.34, la^l decreases, ultimately becomes positive, making 
r.h.s to be positive again. So the jet starts to accelerate at ^ > 6.34 until v vt- 

Now we discuss how various disc parameters and fluid composition of the jet affect 
dynamics of the jet. The jet is being affected by the radiation from the disc, and the radiation 
held above the disc in inhuenced by (3, rhsk, and rukd- The base values of the jet are Vh = 0.014 
and ai) = 0.51 at z^ = 1.5. In Fig. |H1. we show the comparison of u as a function of 2 ; for 
various values of /9 —)■ 0.0 (dashed, black online), 0.5 (solid, black online) and 1.0 (dotted, 
black online). The corresponding terminal speeds {vt) at z = lO'^ with f3 are presented in 
Fig. [HI). As the magnetic pressure increases in the disc i.e., (3 increases, the supersonic part 
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of the jet is accelerated because 4k increases. When magnetic pressure is zero {(3 = 0.0) i.e., 
the jet is thermally and radiatively driven only by pre-shock bremsstrahlung and thermal 
photons, the terminal speed is at around 0.44. And when magnetic pressure is taken to 
be equal to the gas pressure, then vt increases above 0.7. In Fig. [Hb we show the effect of 
accretion rate of SKD on the v prohle and the corresponding terminal speeds are presented 
in Fig. [Hll as a function of rhsk- We see that vt ranges from 0.42 to 0.72 when the rhsk is 
varied from 0.1 to 11.5. The velocity prohle of a thermally driven jet, and a jet driven by 
radiations acted on by rhsk = 1 is similar. Only when the luminosity is close to L-^dd of 
£ —>■ 1, the radiative driving is signihcant. In Fig. [Sb we carry out similar analysis for the 
variation of composition parameter (,^) in the jet, and plot v prohles for jet with ^ = 0.3 
(dashed, magenta online), 0.5 (dotted, magenta online) and 1.0 (solid, black online). With 
the lighter jet v increases, and this is also seen in the vt dependence of ^ in Fig. [Hf. As 
^ increases, high proton fraction makes the jet heavier per unit pair of particles, and the 
optical depth decreases due to the decrease in total number of leptons. So the net radiative 
momentum deposited on to the jet per unit volume decreases, in addition the inertia also 
increases. This makes the jets with higher ^ to be slower. 

For various values of rhsk, Xg changes and therefore not only 4k changes but £ps changes 
too. Infact since Xg is the inner edge of the KD, 4d will change even though rhkd is kept 
constant. In Fig. [9b, we plot vt with the total luminosity i (= 4s + ^sk + 4d), by tuning rhsk- 
As the luminosity of the disc increases the terminal speed increases from moderate values of 
0.44 to high speeds of ~ 0.8 when the disc luminosity is closer to Eddington limit. However, 
4d has limited role in determining vt as has been shown in Fig. [ 9 ) 3 . 

In Fig. [Hthe radiative moment around a super-massive BH is shown to be signihcantly 
higher than that around a stellar mass BH even if the accretion rates (in units of Medd) are 
same. In order to study the effect of the mass of the central object, in Fig. [TOh . we compare 
the V prohle of the jet around IOMq BH (solid, black online) with that around 10 ®Mq BH 
(dotted, blue online). The jets are launched with the same base values {zb = 1.5, Vb = 0.014, 
and Qb = 0.51), and the accretion rates in units of M^dd are exactly same. Although the 
radiative moments around a super-massive BH are signihcantly diherent, yet the v prohles 
diher by moderate amount. To ascertain the cause we plot 3? or relative contribution of 
radiative moments for both the jets in Fig. [TOb . 3? is quite similar for both the BHs close 
to the horizon, but in the range 4 < z < 10 3? around stellar mass BH is higher than that 
around super massive BH. In Fig. [TOb . we compare the Lorentz factor 7 of a jet around 
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Figure 9. Dependence of terminal speeds vt on (a) i and on (b) rhkd for the same base values of jet = 1.5, = 0.014. 

afj = 0.51). The value of rhkd = 1 for (a), and rhsk = 5 for (b). While 0 = 0.5 for both the cases. 


IOMq BH (solid, black online) with a jet IO^Mq (dotted, blue online), launched with hot 
base {zb = 1.5, Vb = 0.19, ab = 0.576) and acted by high accretion rate rhsk = 12. The initial 
jet 7 (= v) is almost same for both the jets, however, due to larger around a stellar mass 
BH, the jet around it slower compared to that around super-massive BH. In this case, the 
terminal Lorentz factor is signihcantly larger for a jet around IO^Mq BH. 

It is clear that jets around stellar mass BH are slower, and lighter jets are faster, but 
what is the maximum terminal velocity possible? We choose to launch jet with maximum 
possible sound speed at the base and very high accretion rate. In Fig. [TTl we plot the 
maximum terminal Lorentz factor or yrmax possible as a function of ^ for IOMq BH (solid, 
black online) and 10®Mq (dotted, blue online), when accretion parameters are rhgk = 12, 
(3 = 0.5 and rhkd = 1- For jet composition ^ = 1.0, the maximum possible terminal Lorentz 
factor 7Tmax ~ 9 for a jet around IOMq BH, but for 10®Mq BH, 7Tmax ~ H. However, for 
lighter jet around stellar mass BH 7rmax ~ 10, but for super-massive BH, light jets yields 
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Figure 10. Comparison of (a) v profile for jets around IOMq (solid, black online) and IO^Mq black holes (dotted, blue 
online). The SKD accretion rate is rhgk = 10. (b) 5R = (E -\-V') jlT with 2 : from an accretion disc around IOMq BH (solid, black 
online) and IO^Mq BH (dotted, blue online). The jet base values for (a) and (b) are 25 = 1.5, Uf, = 0.014, and = 0.51. (c) 
7 profile of a jet around IOMq BH (solid, black online) and 10^ Mq (dotted, blue online), (d) 5R = (£^ + V')I‘1T with 2 . The 
SKD parameter is rhsk = 12. Other disc parameters are rhkd = T /^ = 0.5. And the jet base values for (c) and (d) are Zfj = 1.5, 
Vf, = 0.19, and a^, = 0.576. 


7Tmax ~ few X10. So for light jets ultra-relativistic jets around super-massive BH is possible 
if it is driven by the radiation from a luminous disc. 


5 DISCUSSION AND CONCLUDING REMARKS 

In this paper, we have investigated the interaction of a relativistic fluid-jet with the radi¬ 
ation held of the underlying accretion disc. The accretion disc plays an auxiliary role, in 
other words, the jet-disc connection has been ignored in the present study. In principle the 
accretion rates, and the outer boundary condition, determines the jet states, as well as, 
the radiation held around it. We reduced the number of free parameters to determine the 
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Figure 11. Variation of 7 Tmax with around IOMq BH (solid, black online) and 10 ®Mq BH (dotted, blue online). Accretion 
parameters are rhgk = 12, /3 = 0.5 and = 1. and jet base values are 25 = 1.5, & = 0.5766. 


radiation field, by estimating the temperature, density and the velocity prohle of the var¬ 
ious components of the accretion disc and considered them as generic (Appendix A). The 
radiative intensities of the SKD is determined from the estimated temperature, density and 
the velocity prohle by supplying the accretion rates of rhsk (Appendix |^. The same for KD 
is obtai ned only from the rhkH- We used some known accretion disc solutions from the lit¬ 
erature (IKumar fc Chattopadhva.vll2ni4J ). to obtain a simplifying analytic relation between 
the shock location and the SKD accretion rate (Appendix B). This enables us to reduce 
one free parameter (he. , x.d- Now we use those same accretion disc solutions in the general 
radiative transfer code flMandal fc Chakrabartil 1200811 . and compute the spectra, as well as, 
ratio of the luminosities from the PSD and the pre-shock disc as a function of shock location 
(Appendix B). In other words, the accretion rates determine the accretion solutions and the 
shock location, which determines the disc luminosities, intensities and therefore the entire 
radiation field around it. And then we studied the propagation of a transonic and relativistic 
jet through this radiation held. 
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We noticed that proper relativistic transformation of the radiative intensities from the 
local disc frame to the observer frame is very important and these transformations modify the 
magnitnde, as well as, the distribntion of the moments aronnd a compact object. The PSD 
and SKD are the major contribntors in the net radiative moments and KD contribution is 
much lower than either one of the former two. However, the contribution of PSD compared to 
that of SKD in the net moment, will depend on where the shock forms which is also dictated 
by the accretion rate. One of the interesting fact about the moments due to various disc 
components is that they peak at different positions away from the disc plane. Therefore the 
jet which is initially thermally driven, but is further accelerated by the radiative moments 
from various disc components further down stream. The jets with normal conditions at the 
base, produces mildly relativistic jets with terminal velocities vt ~ few xO.l (Fig. [6]). The 
scattering regime maintains the isentropic nature of the jet, and because we considered a 
realistic and relativistic gas equation of state, the adiabatic index changes along the jet. 
However, radiation not only accelerates but also decelerates if n > Ueq- Although close to the 
jet base (zb) the velocity is low and the radiation held should accelerate, but being hot the 
effect of radiation is not signihcant in the subsonic branch because of the presence of inverse 
of enthalpy term in the radiative term a,, of the equation of motion (equation [15]) . Therefore, 
in the subsonic domain the jet is accelerated as a result of competition between thermal and 
the gravity term. In the supersonic domain the gravity term is smaller than both the thermal 
and the radiative term. If a,. < 0 and la^l ~ at, then the jet can be decelerated (Figs. [7tL,b). 
As the magnetic pressure is increased, the synchrotron radiation from SKD increases and 
it jacks up the flow velocity in the supersonic regime. Increasing the rhgk increases both 
the synchrotron, as well as, bremsstrahlung photons from the SKD, which makes the SKD 
contribution to the net radiative moment more dominant, and therefore increases the v in 
the supersonic part of the flow. The vt increases with (3 but tends to taper off as —)■ 1, 
however, vt tends to increase with rhsk and shows no tendency to taper off. The jet tends 
to get faster with the decrease in protons {i.e., decrease of ^). We do not extend our study 
to electron arid pos itron or = 0 jet, since a purely electron-positron jet is highly unlikely 
(iKumar et a/.ll2013l) . although pair dominated jet {i.e., 0 < ^ < 1) is dehnitely possible. So 
we extend our study to 0.27 ^ ^ 1 jets (Fig. [8^-f). Although the terminal velocity vt 

increases with the total luminosity i, and approaches relativistic values as the disc luminosity 
approaches the Eddington limit. However, the KD plays a limited role in accelerating jets 
(Figs. [9]a,b). Since the radiation held around a super-massive BH and stellar mass BH is 
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different, we compared jets starting with the same base valnes, same accretion rates (in nnits 
of MEdd) bnt aronnd stellar mass and snper-massive BH. If accretion rates are moderately 
high, then jets aronnd super-massive BHs are slightly faster than those around stellar mass 
BH (Figs. [TUh . b). However, if accretion accretion rates are high such that the moments 
from the central region of the disc differ significantly, then ultra-relativistic jets around 
super massive BHs can be obtained compared to stellar mass BH (Figs. ITOb . d). For high 
accretion rates, it has been shown that around stellar mass BH, moments due to PSD and 
the pre-shock disc are comparable, which makes 3? > 1 for a larger distance away from the 
BH (Fig. 0]). This limits the maximum velocity a jet may achieve around a stellar mass 
BH, compared to that around a super-massive BH. This brings us to the issue of maximum 
possible terminal velocity of the jet. We launched jets with relativistic sound speed in order 
to get maximum thermal driving, while increase the accretion rates to very high values. The 
comparison of jet 7 Tmax around a stellar mass BH and that around super massive one as a 
function of shows that jet around super massive BH can be accelerated to 7 t ~ 11 even 
for fluid composition e~ — while that around stellar mass BH 7 ^ ~ few. However, lighter 
jets around IO^Mq BH can be accelerated to truly ultra-relativistic speeds, compared to jets 
around stellar mass BH. 

It is interesting that the radiative moments from different disc components maximize 
at different distances from the jet base, which opens up the possibility of multi stage ac¬ 
celeration of the jet. However, close to the base, thermal acceleration is the main driving 
force that makes the flow supersonic, and thereby negating the effect of gravity. Thereafter, 
the radiative driving accelerates jet to mildly relativistic terminal velocities if the radiation 
held is mild, but intense radiation held can accelerate jets to ultra-relativistic velocities. 
Hence, the thermal driving of the jet at the base, along with the radiative driving further 
out, comprises a complete multi-stage acceleration rather than just the radiation held it¬ 
self. Although the maximum terminal speeds achieved is truly ultra-relativistic, especially 
around super-massive BHs, but at the jet base, which is hot, the disc photons may take 
away energy from the hot electrons instead of imparting its momentum onto it. The jet base 
may also radiate via other processes. These radiations would, in the observer frame, actually 
be howing along the jet and therefore may interact and deposit momentum on to the jet 
further down stream. The jet may also gain energy via free-free absorption. These aspects 
have been ignored and therefore the conclusions might slightly differ. However, it is clear the 
jet can be accelerated to relativistic terminal speed. We would like to simulate the present 
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work of radiation driven relativistic jets aroun d compact objects, similar to our effort with 
galactic outflows (ICliattopadhvav et a/.N2012l) . One of the major difference expected is, if 
the radiative states of the accretion disc changes, how does it affect the jet and what are the 
timescales in which these changes are expected to be observed in the jets. These advances 
of our work are underway and will be reported elsewhere. 
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Figure Al. Velocity and dimensionless temperature ©a estimated for accreting flow for angular momentum A = 1.7. Here 
ro = 3500 and = 0.01. 
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APPENDIX A: ESTIMATION OF TEMPERATURE AND VELOCITY OF 
SKD AND PSD 

In the present analysis, the focus rests upon outflows and the accretion disc plays an auxiliary 
role, as radiation from the various accretion disc components affect the jet dynamics. Since 
the accretion disc emission depends on the flow variables of the disc, we need to know the 
temperature and velocity distribution. The method to estimate the velocity distribution of 
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SKD and PSD are given below. We estimate the velocity from geodesic eqnations, which 
shows that the covariant time component of fonr velocity i.e., Ut is a constant of motion. 
From = — 1, hnd Now, 

—UtU^ = = 7 ^ = 7 ^ 7 |=sqnare of Lorentz factor, 


where. 


= 


1-^9?’ 




v^„ = 


—Ufj)U'^ (r — 1) 


V 


2 1 ' <j>a. 
(pa 


UtW- 


(Al) 


Here, is the radial 3-velocity measnred by a co-rotating observer, while the radial 3- 
velocity is Vra, and are dehned as 


v^ = 

ra 


-UrU 


UtW- 


= iVr: 


The snffix ‘a’ stands for either PSD, SKD or KD. For KD v^kd = "dkd = 0 and v^kd is fhe 
Keplerian azimnthal velocity. We chose onter bonndary conditions (da = doa, and Aq at 
r = ro), this relation allows ns to calcnlate Ut for SKD, 


MX (i )^_ 


X'4 - X - 1)A§' 
And then from eqnation IA2l we obtain d(r)sk, 


(A2) 


dsk = 


1 - 


(r — l)r^ 


1/2 


(A3) 


X- [{r -l)XXuX. 

Here, the eqnations are expressed in geometric nnits 2G = Mb = c = 1. One mnst note from 
eqnation flA3p that as r ^ 1, d —)■ 1, althongh SKD in presence of shock, does not extend 
npto the horizon. However, one shonld also remember since the velocity is obtained from 
geodesic eqnations, it is slightly over estimated becanse the pressnre gradient terms were 
ignored while estimating the velocity. This would under estimate the radiative moments 
slightly and hence our results are believable since there is no over estimation of jet and 
radiation interaction. To get the velocity distribution for the PSD, we assume the shock 
compression ratio is 3, so the ratio of post shock (suffix -|-) and preshock (with suffix —) 
velocities are d+ = d_/3 at r = x^, with this boundary condition we recalculate the constant 
of motion and then use it to estimate the d distribution for PSD, 

,2 1 1/2 


d = 

^ps 


(r - l)r" 


{-■“-|(r-l)V]KuJ ■ 

In Fig. flAlh ) we have plotted the velocity distribution of SKD and PSD, where the suffix 
‘a’ signihes either sk or ps. 
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A1 Density and temperature of SKD 

The SKD accretion rate equation is 

Msk ~ (AS) 


here is the local height of SKD. If the outer boundary of the accretion disc is at tq, the 
height at outer boundary Hq, the radial four velocity Uq, the dimensionless temperature ©o 
and the density po, then we have 


Psk _ 

Po Msk^-^sk 


(A6) 


Here, u^y- is obtained from equation IA3I and Hgy = rcotOgk + do. 

Assuming fixed T and moderate radiative loss, then the dimensionless temperature or 
©sk = Psk/PskC^ is given by. 


^ = Q ( 


©sk — ©0 


r -1 


Po 


V <krH J 


(A7) 


Equations flA7l and IA6I) are used to compute the intrinsic synchrotron and bremsstrahlung 
intensity of the SKD. Using the expression of = {1 — equation IA71 the 

SKD and PSD temperature distribution can be obtained. In Fig. flAlb ). the dimensionless 
temperature distribution ©a is plotted for both SKD and PSD. 


APPENDIX B: RELATION BETWEEN M,y AND X, 


In Fig. (iBlla) red triangles represent the locations of shock for a given rhsk- Other parameters 
are ro = 3686, Pq = 1-928 x 10“^, ©o = 9.811 x 10“^, Aq = 1.7 ar id the viscosity parameter 
a = 0 .001. In the figure, Xg is computed using the methods of iKumar fc Chattopadhvav 
(l2ni4J ) for the gi ven values of a, rhsk, rp, po, Qn and Aq. We fitted a curve using these 
data generated bv iKnmar fc Chattopadhvayl (120141 ) and expressed Xg as a function rhsk, the 
explicit form of which is 


Xg = 64.8735 - 14.1476msk + 1.24286mi - 0.039467m: 


^sk* 


(Bl) 


Therefore, while calculating the radiative moments, Xg is no more a free parameter, but is 
estimated using equation (IBID . 
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Figure Bl. (a) Shock location Xs as a function of rhsk (filled triangles, red online). The fitted function given by equation m 
(solid curve, black online). Parameters fixed to get the data (filled triangles) are ro = 3686, i9o = 1.928x 10 ^, ©o = 9.811 x 10 ^, 
Aq = 1.7 and a = 0.001. (b) x a function of Xg for IOMq (filled square, black online) and IO^Mq BH (filled circle, blue 
online). The corresponding fitted functions for IOMq BH (dashed, black online) and 10^ Mq BH (solid, blue online) presented 
in equations |[B3l[B4ll . (c) Spectral index as a function of rhgk plotted for the same data points of Fig. a. for IOMq BH (filled 
square, black online) and 10^ Mq BH (filled circle, blue online). 


Bl Obtaining relation between shock location and ratio of pre-shock and 
post-shock Inminosities 

V driven jets were stuc 


In previous works in which radiative 


2nnnH 


2002 


Chattopadhvay et al. 


20041 : 


ied fjCliattopadhyav fc Chakrabarti 

Chattopadhvayl 200^ . the luminosity from PSD 


(£ps) was supplied as a free parameter, but in the present work, we calculate this from 
spect ral modeling. KD produces the well known thermal radiation fjShakura fc Snnvaev 


19731), while SKD emits via bremsstrahlung and synchrotron emission. For the same set 


of outer boundary conditions, i.e., a, rhsk, ro, ho, 0o an d An and rhwH at rp, we so 


general radiative trans f er equ ations similar to those by 


Chakrabarti fc Titarchuk 


ve the 


199.4 : 


Mandal fc Chakrabarti! (120081 ). For the disc parameters of Fig. [B]la, the shock location var¬ 


ied between 51.25 ^ Xg ^ 12.6 as on e varied 1.1 ^ rhsk ^ 7, k e eping all other boundary 


conditions same. The flow solutions of 


Kumar fc Chattopadhvayl (120141) were obtained in a 
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unit system where Mb = 1, so the actual luminosity and spectra will depend on the mass 
of the central object. PSD being optically slim will inverse-Comptonize radiations coming 
from pre-shock disc (SKD and KD). Let the luminosity from SKD and KD be denoted as 
ips and f'kd- Then the ratio of luminosities is dehned as y = ipg/ (4k + 4d)- In Fig. (iBl lb) we 
plot y with Xs for the same set of data points as in Fig. fiBlla) and considering rhkd = 1, we 
obtain the spectra and luminosities of each disc component for Mb = IOMq (black, squares) 
and Mb = 10®Mq (blue, circles). 

The behaviour of y for 10 ®Mq is different from that of IOMq because synchrotron cool¬ 
ing is very efficient for a stellar mass black hole. For a IOMq black hole, the post-shock 
luminosity increases initially as the post-shock flow is still very hot but as shock location 
moves close to black hole, rhsk is very high and pre-shock synchrotron cooling becomes large 
enough to reduce the ratio (y). On the other hand y keeps on increasing for 10®Mq because 
Comptonization due to Keplerian soft photons is the most efficient cooling process and as 
shock moves in, the supply of hot electron increases (as rhsk increases) which enhances the 
post-shock luminosity. In Fig. flBifcl we plot the variation of spectral index with rhsk- We see 
that as rhsk increases the spectral states becomes harder because supply of hot the electron 
increases. Moreover, the spectral index for a super-massive black hole is generally softer than 
a stellar mass black hole as the PSD of stellar mass black hole is relatively hotter. We htted 
the plots of Fig. fiBllb) with analytic functions (the constants are written correct up to three 
decimal points) given by 

yg = 25.944 - 1.667a:, + 3.992 x lO'^a;^ - 3.199 x lO'^a:^ (a:, > 12); (B2) 

yg = 1.449 + 2.336a:, - 0.127a:2. (a:, 4 12) 

yi = -2.525 + 0.913a:, - 4.438 x lO’^a:^ + 6.522 x lO’S^ (x, 4 29); (B3) 

yi = -2.914 - 1.622 x lO'^x, + 7.265 x + 1.278 x 10"^x^ (x, > 29). 

Here, yg signihes the ratio of PSD to SKD-I-KD luminosity for Mb = 10®Mq (blue, solid 
curve Fig. [B]lb), while xi is the ratio of luminosities for Mb = IOMq (black, dashed curve 
Fig. (Bjlb). Therefore we supply rhgk and parameters at the outer boundary and estimate 
the hsk, ©sk and psk- up to the x, (obtained via equation IBlIi and these solutions we obtain 
4k- The 4d can be estimated from rhkd, ro (outer edge) and x, (inner edge). So from the 
preshock luminosity (4k + 4d) we can estimate the 4s using either expression flB4p or flB3jl 
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+ log(m^cVk) 

*log(mpCVk) 

Figure Cl. (a) Comparing adiabatic index F as a function of temperature T, for RP EoS (solid, black online and long dashed, 
magenta online). CR EoS is presented by dashed curve (red online) and dotted (blue online) curves. The two separate cases 
^ = 0 and ^ = 1 are marked on the curves. Temperature corresponding to electron and proton rest mass are also marked on 
the T axis, (b) T-p as a function of © for e~ — p"*" flow, but for RP EoS (solid, black online) and CR EoS (dashed, red online). 


depending on the central mass. In this work IOMq is considered as a representative of stellar 
mass BH and IO^Mq is considered as super massive BH. 


APPENDIX C: EQUATION OF STATE 

The EoS for relativistic fluid is obtained by integrating the relativistic energy of fluid par¬ 
ticles followi ng a relativistic Maxwe ll -Boltzman di st ribution in th e morn entum space as was 
obtained by Chandrasekhar (1938); Synge ( 1957 1: Cox fc GinlJ ( 1968 1. For single species 


fluid the different forms of the EoS obtained by the above authors are as below, 




P; ecG = pc 


K2{1/Q)J 


,(C1) 


4A2(1/0) ’ " ^2(1/0) 

where, e represents local energy density of the flow, 0 the measure of temperature and 
suffix C, S, and CJ signihes Chandrasekhar, Synge and Cox & Giuli, respectively. The K’s 
are modihed Bessel’s functions of second kind and respective indices indicate their degree. 
We recall the recurrence relation Km+i{x) = Km-i{x) -|- 2mKm/x, and obtain 
e+p 


he, = 


pc^ 




(C2) 


Therefore, all the forms of exact EoS for relativistic gas presented above, are equivalent, and 
let us denote the EoS represented by equation (ICIH as ‘relativistically perfect’ or RP EoS. 
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Figure Dl. Estimated width of the jet by (crosses) by matching the thermal, gravity and radiative forces along r, for jet 
base values 25 = 1.5, = 0.50544, V}, = 0.00345 and disk parameters rfisk = 7 and rhkd = 1- The conical cross-section 

approximation is shown by a fitted line (dashed, red online). The disc top surface and disc components (PSD and SKD) shown, 
is estimated from accretion rates supplied. 


The multispecies approximate EoS used in this paper, may be called Chattopadhyay-Ryu 
or CR EoS (equation [TU]), not only mimics equ ation fICip very well, but also satishes the 


fundamental inequality obtained by 


Taub 


fjl948[l . Taub showed from Erst principle that any 


EoS for dilute relativistic gas have to satisfy a fundamental in-equality given by 


Tt = 




^ 1 . 


(C3) 


By following equation flT^ . we can calculate T for any EoS. In Fig. (Of a) we plot adiabatic 
index T as a function of temperature T, for e~ — p"*" (marked ^ = 1) and e~ — e"*" (marked 
.^ = 0) flow. For e~ — flow, the solid curve (black online) is the T with RP EoS, and 
dashed curve (red online) is the T with CR EoS. And for e~ — flow, T with RP EoS is the 
long dashed (magenta online) curve, and that due to CR is dotted (red online) curve. In Fig. 
fICllb) we compare the Taub function as a function of 0 for e~ — flow described by 
RP EoS (solid, black online) and CR EoS (dashed, blue online). Both the EoS comfortably 
satisfy the in-equality. 
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APPENDIX D: ON THE JET GEOMETRY 


In this paper we assumed the flow surface of the jet to be conical having a small opening angle 
and maintains it throughout the spatial extent. Here we test how good is the approximation. 
To approximately locate the lateral extent of the jet, we balance the pressure gradient term, 
the gravity term and the radiative term along r direction. We assume that dp/{pdr) at the 
jet edge is equal to dp/{pdz) on the axis. So 
IdP 


2re 2 ^ 

V dz z, 


ap. = —TP ~-1 —T- + - I (Dl) 

p or T \ V dz z) 

Similarly, approximating k. k. 0, the component of the radiative term along r can be 
evaluated as (see, equation 3a of Chattopadhyay 2005), 

-VV^^) (D2) 

The gravity term is 

+ (D3) 

In Fig. (IDL) we plot the location where |ap^| = \ar^ + a^^l for a radiation held due to disc 
parameters rhgk = 7, rhkd = 1 and jet characterized by zi, = 1.5, = 0.50544, vi, = 0.00345. 

The accretion disc upper surface is also shown in the figure. The locations where |ap^| = 
l^r-r P ^gr\i is shown by crosses. The approximated boundary of the jet (dashed, red online), 
shows consideration of conical jet flow geometry is a fairly good assumption. Although, this 
assumption is reasonable only when the jet opening angle is small. This is because for large 
opening angle, jet material at the edges would be spun up by which may contribute in 
spreading the jet. 
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